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Low temperature line-width broadening in optical-conductivity spectra of the 
off-center rattling phonons in type-I clathrate BagGaigSuso 
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With a teraiiertz time-domain spectrometer (0.3 - 3.0 THz) we have measured the optical conduc- 
tivity of the type-I clathrate BagGaieSnso at temperatures from 300 K down to 7 K. Independent 
six spectra superimposed on the Drude conductivity are identified to infrared active vibrational 
modes of guest Ba ions and the cages. While the spectra of five higher-frequency modes depend 
hardly on temperature, the lowest-lying spectrum with a peak at 0.72 THz due to the Ba(2) ion's 
off-centering vibration in the oversized cage changes with temperature characteristically. With low- 
ering temperature, the spectral shape of this so-called rattling phonon continues to become so broad 
that the line-width amounts to be comparable to the peak frequency. Furthermore, below about 100 
K, the single broad peak tends to split into two subpeaks. While this splitting can be explained by 
assuming a multi-well anharmonic potential, the strong enhancement of the line-width broadening 
toward low temperature, cannot be understood, since the Boltzmann factor generally sharpens the 
low-temperature spectra. 
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Recently novel phenomena related with local anhar- 
monic phonons observed in cage-like materials such as 
clathrates[l|, y] have attracted much interest. In partic- 
ular, in relation to their potential applications to ther- 
moelectric devices, it has been studied extensively how 
the so-called rattling phonons coined for the local anhar- 
monic phonons of guest ions in polyhedral cage could be 
responsible for such electric and thermal conductivities 
and the specific heat. 

In order to clarify the dynamical properties of these 
phonons, we reported in the previous paper[3| the 
optical-conductivity spectra in the terahertz frequency 
region for the type-I clathrate BagGaieGeao (BGG), i.e., 
a p-type semiconductor with a low carrier-density. It 
has turned out that the vibrations of Ba^+(2) ions in the 
tetrakaidecahedra framework are quasi on-centering in an 
anharmonic potential, indicating characteristic features 
as the softening of the peak frequency and the sharpen- 
ing of the line- width toward low temperaturel3|,|j]. These 
two features of the on-center rattlings are explained by a 
local anharmonic potential (AP) model. An on-center 
type anharmonic potential has a non-equivalent level 
spacing contrary to the harmonic potential; the level 
spacing becomes larger as the level becomes higher. As 
the effect of the Boltzmann factor, many and higher levels 
participate to optical transitions at higher temperature, 
so that the peak frequency becomes higher and the to- 
tal line-width becomes broader. On the other hand, at 
low temperature, local phonons are almost condensed at 
the ground level and hence the relevant optical excitation 
energy becomes lower accompanied with a suppression of 
the line-widthQ. 



As an extreme case of anharmonicity in the local po- 
tential, one can consider an off-centered anharmonic po- 
tential. In general such a potential has multiple lo- 
cal minima and may cause the low-lying rattling vi- 
brations with a large displacement. Through couplings 
with those modes, one expects strong effects on elec- 
trons and acoustic phonons. In fact, the type-I clathrate 
BGS, which is isostructuralpj, |5[ to BGG, is found to 
show the off-centered displacement of Ba^+(2) by the 
single-crystal x-ray diffraction measurement [ij. More- 
over there have been revealed anomalous properties such 
as a glass-like behavior in the specific heat and thermal 
conductivity m |2| as well, and a resistivity niinimum[2|. 
Those behaviors have been discussed in relation to the 
off-centered rattling phonons studied by Raman scatter- 
ing measurement [6| . Infrared-active optical conductivity 
measurements, complementary to the Raman scattering, 
are expected to shed more light into the dynamical be- 
haviors. Theoretically, in Rcf. [J] for example, it is pre- 
dicted that the off-centered anharmonic potential would 
produce a line splitting in the optical conductivity at low 
temperature together with a softening and/or hardening 
of the peak frequency. 

In this letter, we measure the optical conductivity of 
infrared-active optical modes of type-I BGS by the ter- 
ahertz time-domain spectroscopy (THz-TDS) to directly 
identify the off-centering rattling phonon and reveal its 
dynamical properties. We discuss the observed results 
based on the AP model [J|. The observed spectra, which 
show the line-splitting just as expected from the model, 
reveal an unexpected, drastic line-width broadening to 
occur towards low temperature, which is quite opposite 



from the usual thermal behavior. 

Single crystals of the type-I BGS doped n-type carriers 
are grown by a self- flux method [l|,l4|. To obtain sufficient 
transmitting signals through the conductive sample, sin- 
gle crystal disks of 3 - 5 mm in diameter are polished down 
to about 13/xm in thickness by using diamond abrasive 
films. The present measurements covering a frequency 
range of 0.3 - 3.0 THz are carried out with the commer- 
cial spectrometer (RT-20000, Tochigi Nikon Co. Ltd) 
which uses a standard technique for the transmission con- 
figuration and the temperature range from 7 K to room 
temperature [si [D] . 

Figure [1] shows the real part, 0-1(0;), of the complex 
conductivity for data at temperatures from room tem- 
perature down to the lowest 7K. Just as previously 
found in the isostructural BGGl3[, the data demonstrate 
four broad peaks superimposed on almost frequency- 
independent Drude-like contributions (see horizontal dot- 
ted lines in (b) and (c)) caused by doped carriers. The 
latter contribution monotonically increases with decreas- 
ing temperature. This metallic behavior is contrary to 
the semiconducting behavior of p-type BGG in terahertz 
frequency region [3J. 

The observed four peaks are assigned as six indepen- 
dent modes {vi — vq) by the first-principle calculations 
for the infrared active optical modes of Tiu symmetry [6|. 
To reproduce each spectral shape, we define the fitting 
curve as 



aijit{uj) = 0-Q + 



{2/tt)S,uj^T, 



{c^l.-Lo^r + icor^y 



(1) 



where the first term ctq is the dc conductivity (w = 
0) of the carrier contribution, and the second term is 
a Lorentz-type conductivity for the phonon mode Vi. 
There, Si, coqa, and Ti are a spectral weight, a resonant 
angular frequency and a relaxation rate (or a damping 
factor), respectively. We apply this equation to fit all the 
data and demonstrate the fitting results at room temper- 
ature and 7K in Fig. [TJb) and (c). Fitting parameters 
are listed in Table HI 

Firstly, the lowest peak labeled vi is assigned as the 
mode of the Ba^"'"(2) in the oversized tetrakaidecahedra 
cage. In order to reproduce the line shapes, this peak 
is fitted by superposing two Lorentzian curves, Lia and 
Lib. The peak frequency of wq,! in Table U is defined 
as the peak position of the superposed Lorentzian curves 
(see, the arrows in Fig. [l](c)). The relaxation rate Fi of 
1^1 is evaluated from the difference between the frequen- 
cies of half maximum value for the curves. The vi mode 
corresponds to the vibration within the plane perpen- 
dicular to the fourfold inversion axis, i.e., the off-center 
rattling phonon mode in the oversized cage. 

Secondly, the small peak around 1.3 THz labeled ^^2 is 
also assigned as the mode of the Ba^"'"(2) in the oversized 
tetrakaidecahedra cage. This mode is the vibration along 
the fourfold inversion axis. 
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FIG. 1: (color online) (a) Temperature-dependent real part 
of complex conductivity spectra (Ti(a;) of type-I BGS. (Note 
the horizontal axis is frequency (= uj/2n) in THz). (b) and 
(c) Experimental data (open circles) and fitting results (lines) 
at room temperature and 7K. Details are described in text. 
For convenience, ITHz = 33.3cm"^= 48 K= 4.14 meV. 



Thirdly, the third peak spectrum around 1.7 THz is re- 
produced by superposing two Lorentzian curves, where 
the larger one is assigned as the mode labeled 1/3 of 
Ba^+(1) in the smaller dodecahedra cage and the smaller 
one as the mode labeled 1^4 of the cage. 

Finally, the peak around 2.2 THz is also fitted by su- 
perposition of two Lorentzian curves, and the both curves 
are assigned as the cage modes labeled 1/5 and i/q. It is 
to note that the spectral weight Si for all the modes are 



TABLE I: Fitting parameters at 296 K (7K) of type-I BGS. 
Tiie guest ions in the smaller (2a) and oversized (6d) cages 
are defined as Ba(f) and (2), respectively. The definitions of 
ojo.i and Fi are described in text. 
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independent on temperature within an experimental er- 
ror. 

As can be seen in Fig. [T] (a), all the spectral shapes 
(i/2 — i^e) except the lowest i^i do depend little on tem- 
perature. For the peak frequencies, the 1/2 mode soft- 
ens by 1.5% and the 1/3 mode also softens very slightly 
as temperature decreases, while other higher modes of 
cages become constant or hardened. These behaviors are 
essentially similar to those in BGG[3[. 

Hereafter we will concentrate on the i^i mode of the 
rattling phonon that is drastically temperature depen- 
dent. As seen from Fig. [T] (a) and also shown by circles 
in Fig. [2] (a), the peak frequency of ui mode ojo.i slightly 
softens by 7 % toward 100 K and turns to harden at fur- 
ther lower temperature. This tendency has never been 
observed in on-center type rattling phonons in BGG|3[. 
The line-width Fi indicated by marked diamonds in Fig. 
[2] (b) increases almost by a factor of 2 from 0.3 THz at 
room temperature to 0.6 THz at 7K. 

For clarity. Fig. [3] expands the spectral shapes at room 
temperature (diamonds) and 7 K (circles) already shown 
in Fig. [TJ At 7 K, the i^i spectrum clearly deforms 
around 0.6 THz, and has been well reproduced by as- 
suming the double peaks at about 0.5 and 0.7 THz. 

Based on the analysis using ID-AP model, the splitting 
at low temperature is well understood using the concrete 
vibrational levels in the double- well potential V{x) = 
l/2kx^ + l/4Aa;'* shown in Fig. U) The parameters are 
fixed to reproduce the frequencies, wo,iA=0.52THz and 
wo,iB=0.74THz for the double peaks. Using the value 
of the Ba mass, potential parameters are uniquely deter- 
mined; k = -2.9(kgs-2) and A = 3.7 x 102i(m-2kgs-2). 
The features of the calculated levels are as follows. 
Firstly, the n = and 1 levels of are almost degenerated 
with a tiny energy difference hujio = 0.059 meV. Secondly, 
the n = 2 and 3 levels are well separated due to the mod- 
erately deep potential well. From the ground level n = 0, 
hu}2o = 2.17 meV (=0.52 THz) and htoso = 3.05 meV 
(=0.74 THz), respectively. These characteristic low- lying 
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FIG. 2: (color online) Temperature dependence of (a) the 
peak frequency wo,! (circles) and (b) total line width Fi (dia- 
monds) for vi mode of type-I BGS. The resonant frequencies 
(i^o.iA (open triangles) and ojo.ib (open squares)) and the re- 
laxation rates (Fia (open inverted triangles) and Fib (open 
pentagons)) are evaluated from the fitting for vi mode. 



levels from ?i = to 3 can produce double peaks at 
low temperature as follows. At sufficient low tempera- 
tures, phonons are condensed to almost degenerate levels 
of n = and 1. Therefore the optical transitions occur 
mainly from n = to 3 and from n = 1 to 2 due to a par- 
ity requirement of the wave function. These transitions 
result in forming the double peaks with the excitation 
energies, U!2i ~ ^^^20 = 0.52 THz and 0^30 = 0.74 THz. 
With these excitations and the Lorentzian curve (eq. 1), 
the data are best fitted as shown in FiglJl leading to 
FiA = 0.46 THz and Fib = 0.32 THz just, respectively 
corresponding to F30 and T2i. 

However, the significant broadening observed in the 
rattling phonon spectra toward low temperature cannot 
be reproduced by our calculations based on the ID-AP 
model or its extension to more realistic potentials of 
higher order and dimension. (The details will be pub- 
lished elsewhere.) 

In conclusion, the temperature dependence of the 
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FIG. 3: (color online) Expanded views for the spectral shapes 
of vi (in Fig[T]) at room temperature (diamonds) and 7K 
(circles) and fitting curves at 7K (lines). Experimental data 
are analyzed by subtracting the fitting results for neighboring 
phonon modes {1^2, 1^3) and the flat Drude conductivity from 
the raw data. Thin lines are the Lorentzian curves of Lia and 
Lib) with the relaxation rates Fia and Fib- Thick line with 
Fi,7K is obtained by superposing these two curves. 



conductivity spectra has been determined for the off- 
centering rattling phonon in the oversized cage of type-I 
BGS. Using the ID-AP model, the spectral splitting at 
low temperatures can be explained, qualitatively at least, 
as optical excitations from highly degenerate ground and 
the first excited levels to the second and third ones. 
Then the spectra are well reproduced by superposing two 
Lorentzian shapes. There remains a question, however, 
on the present finding that the line-width, a direct mea- 
sure for the relaxation rates of the rattling mode, in- 
creases continually with decreasing temperature, which 
would be quite opposite to the spectral sharpening usu- 
ally expected from the Boltzmann factor. 

So far we have discussed our spectra with the ID-AP 
and its extended models treating only the anharmonic 
optical modes. All the other possible interactions are 
neglected at all; the dipole-dipole interaction between 
neighboring rattling phononslsl, |9| and the interactions 
with acoustic phonons propagating along the cage net- 
work. The interaction with the conducting charge carri- 
ers has also been discussed [10|. Since the present spec- 
troscopic technique provides us a phonon-mode selective 
window, further experimental and theoretical studies are 
highly encouraged to delineate how the ofF-ccntcr rattlers 
are interacting with these excitations. 
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